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Abstract 

General properties of charge carrier transport in disordered organic materials are dis- 
cussed. Spatial correlation between energies of transport sites determines the form of 
the drift mobility field dependence. Particular kind of spatial correlation in a disor- 
dered material depends on its nature. Mobility field dependences have to be different 
in polar and nonpolar materials. Different methods of mobility calculation from the 
shape of photocurrent transient are analyzed. A widely used method is very sensitive 
to the variation of the shape of the transient and sometimes produces results that 
effectively masquerade the true dependence of the mobility on electric field or trap 
concentration. Arguments in favor of the better, more reliable method are suggested. 
Charge transport in materials containing charged traps is considered without using 
the isolated trap approximation and this leads to qualitatively different results. They 
indicate that the effect of charged traps can hardly be responsible for experimentally 
observed transport properties of disordered organic materials. 
Keywords: charge transport; correlated disorder; computer simulations 

INTRODUCTION 

Any feasible electronic application of a polymer material implies that the material can con- 
duct an electric current via some kind of charge transport mechanism. Velocity of charge 
carriers (electrons or holes) may be an important parameter per se effectively determining 
operating characteristics of the device or it could influence other important parameters such 
as carrier injection or recombination efficiency. Some degree of structural disorder is present 
in most polymers currently used in real life applications, so they naturally belong to the 
general category of disordered organic materials. At the same time disordered low molecu- 
lar weight organic materials often demonstrate similar general features of charge transport; 
for this reason in this article various properties of charge transport in disordered organic 
materials are mostly discussed without specific emphasis on polymeric nature of a material. 

For many years the study of charge carrier transport in disordered organic materials was 
mostly motivated by hope to use them in various electronic and optoelectronic devices. At 
the same time this problem is an interesting and sometimes even intriguing area of modern 
condensed matter physics. Good review of the area may be found in recent monographs.^! 
We consider here only materials which are insulators under usual conditions and attain 
short-lived conductivity after injection of charge carriers under the action of laser pulse or 
electric discharge and completely leave aside such highly conducting materials as doped poly- 
acetylene, polyanyline etc. Typical materials include low molecular weight organic glasses, 
polymers doped with aromatic organic molecules (usually amines or hydrozones), polysi- 
lanes, poly conjugated polymers, and other compounds. This diverse variety demonstrates a 
lot of common transport features unambiguously indicating a common transport mechanism. 

Transport properties of materials are usually characterized by carrier drift mobility \i and 
its dependence on applied electric field E, temperature T, and other relevant parameters. 
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The most direct method of the mobility measurement is the time-of-flight method (TOF) 
where carriers are generated at the vicinity of one electrode of the plane capacitor (organic 
material fills the space between electrodes) and then drift to the opposite electrode under 
the action of applied uniform electric field E. Analyzing the shape of the photocurrent 
transient, one can calculate the transit time tt? (the moment when current begins to drop, 
thus indicating arrival of carriers to the opposite electrode) and then the mobility 

Ai = -«— (1) 
E Vot tr 

where L is thickness of the transport layer, v is the average velocity of charge carriers, and 
Vo is the applied voltage. Dominant majority of experimental data cited in this article was 
obtained by the TOF method. 

Let us briefly summarize major experimental findings offering important clues on the 
nature of charge transport in such materials.^"^ 

1) Molecularly-doped polymers offer a natural possibility to study dependence of the mobil- 
ity on dopant concentration. It was found that 

In^w -2R/Ro, (2) 

where R is the average distance between dopant molecules. This relation suggests that the 
transport occurs as a series of hops between localized states (transport sites) originating, 
in this particular case, at the dopant molecules. However, in some exceptional cases the 
radius Ro of the wave function of a transport site is unreasonably largeP 

2) Mobility increases exponentially with the increase of temperature. There is no consensus 
on the form of this dependence: some authors favor a simple activation dependence 

ln/i oc — 1/T, (3) 

while others argue that the dependence 

ln/i oc -1/T 2 (4) 

is more accurate. It is generally accepted that an exponentially strong temperature 
dependence may be most naturally explained by suggesting that the major factor gov- 
erning charge transport in disordered organic materials is an energetic disorder: random 
locations and orientations of molecules in the bulk of the material produce random fluc- 
tuations of site energies C/j. Then the dependence /x(T) should be connected with the 
density of states P(U). There is no reason for P(U) to maintain the same form in different 
materials and this may be the reason for different forms of /x(T). 

3) The most surprising fact is the mobility field dependence. Indeed, for carrier hopping 
between spatially localized sites we should expect 

\nfioceER/k B T, (5) 
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just estimating a typical drop in the difference of carrier energy between two transport 
sites separated by the typical intersite distance R. Yet instead of eqE]an almost ubiquitous 
Poole-Frenkel (PF) dependence emerges in the materials 



ln/i oc v E, (6) 

leading to much stronger field dependence of the mobility in weak field. Uncertainty in 
the form of dependence /x(T) is mostly originated from the difficulty to measure the quasi- 
equilibrium (nondispersive) mobility in a wide temperature range. In striking contrast 
to this limitation, dependence fl6]) was measured for some materials in really wide field 
range (an unrivaled example is ref. [5] with E varying from 8 x 10 3 to 2 x 10 6 V/cm). 

For a long time an explanation of the PF dependence has been considered as a major problem 
for the charge transport theory. Since the early paper of Gill, who suggested an empirical 
formula 

H = Hq exp 7 Qj _ JL^ (yE - y/W^j , (7) 

based on the study of charge carrier transport in polyvinylcarbazoleP it was very tempting 
to attribute the field dependence of fi to the influence of charged traps. In the case of 
trap-controlled charge transport with charged particles (having charge of the opposite sign 
to carriers) serving as traps, an applied electric field E leads to the decrease of activation 
energy of the carrier escape from a trap. In this case 

(here £ is a dielectric constant) and calculated values of 7 are usually close to the measured 
ones by the order of magnitude, though discrepancies by a factor of 2-3 are typical. Yet this 
classical PF mechanism attained a lot of criticism because there is no evidence for charged 
traps to be a common constituent in very different classes of disordered organic materialsP 
A first really successful transport model was the Gaussian Disorder Model (GDM) Pin the 
GDM transport sites are arranged on a regular lattice and are assigned site energies Ui drawn 
independently from a Gaussian distribution having variance a 2 . This model explained many 
features of charge transport in disordered organic materials, among them the temperature 
dependence of mobility, transition from the quasi-equilibrium transport to the dispersive- 
like with temperature decrease, and others features. Still the explanation of the PF field 
dependence remained an unsolved problem because the GDM can reproduce this dependence 
in the limited field range only, not significantly wider than 3 x 10 5 — 1 x 10 6 V/cm, in 
weaker field dependence fi(E) has the usual form (jSJ). Moreover, the only reason for the 
approximate PF dependence in the GDM is the Miller- Abrahams hopping rateP There is no 
reliable evidence that this particular rate is the universal hopping rate in a wide variety of 
disordered materials. For this reason it was highly desirable to retain all useful features of the 
GDM and yet incorporate some new feature permitting to explain the PF field dependence. 
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CHARGE TRANSPORT IN POLAR MATERIALS 



A history of the GDM modification that eventually lead to the successful explanation of 
the PF dependence ([6]) began with the study of dipolar contribution to the total energetic 
disorder in organic materials. Computer simulation of the distribution of energies for a carrier 
moving through a simple cubic lattice with sites occupied by randomly oriented dipoles (thus, 
the energetic disorder in this case originates from the electrostatic charge-dipole interaction) 
was carried out by Dieckmann, Bassler, and Borsenberger^ they found that P(U) has the 
Gaussian form for high concentration of dipoles and goes to the Lorentzian one for low dipole 
concentration. They calculated the energy variance a 2 = (U 2 ) and found that typical values 
of <7d are 0.05-0.1 eV which means that the dipolar contribution should be a significant part 
of the total a (according to the GDM^cx ^ 0.1 eV is needed to fit experimental temperature 
dependence of /i). Later it was found that for this model of dipolar glass (DG) function 
P(U) could be calculated analytically.^ For a simple cubic lattice 

a d = 2.35^ (9) 

ea z 

(here a is a lattice scale, p is a dipole moment, and c is a fraction of sites occupied by dipoles) 
and the distribution of U, while indeed having the Gaussiam form for c ±2 1, has the long 
tail P(U) oc 1/U 5/2 for c< 1. 




Figure 1. Distribution of site energies Ui in 50 x 50 x 50 sample of the simple cubic lattice 
model of dipolar glass. Black and white spheres represent sites with positive and negative 
values of Ui, correspondingly, while the radius of a sphere is proportional to the absolute 
value of Ui. Sites with small absolute values of U (less than 3ep/ea 2 ) are not shown for the 
sake of clarity. 

The most significant property of the DG model is a strong spatial correlation in the 
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distribution of U (r) and the energy correlation function decays very slowly with distance^ 



C(r) = (U(r)U(0)) « 0.74^-, r > a 



(10) 



(here angular brackets denote statistical average). In a correlated Gaussiam medium the 
conditional probability P(U\Uq) which is the probability density to have energy U at the 
site r if the site energy at the reference site r*o is Uq has the formP^ 



P(U\U 



1 
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Keeping in mind that C(0) = a^, we see that equation ITTl expresses the tendency of sites that 
are spatially close to have close energies, in striking contrast to the GDM case (see FigHJ). 




Figure 2. Computer simulation of the cluster distribution on size (squares - the GDM, 
triangles - DG model with same a). Here N is the number of sites in a cluster, N c is the 
number of clusters and N tot is the total number of sites; cluster is defined as a connected set of 
sites having energies greater than some (arbitrary) threshold energy (3.7a in this particular 
case). 

Another manifestation of the correlated nature of DG is an exact result for the dependence 
of the distribution of average domain energy Uy = (U) v on domain size V. For a Gaussian 
random medium without any spatial correlation the distribution has the Gaussian form with 
variance 

a 2 v = a 2 a s /V, (12) 

while for spherical domains in DGP^ 

air = 12a 2 d a/5R, (13) 
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where R is the radius of the domain. The difference between (fT2l) and ({TBI suggests that 
large domains are much more common in dipolar matrix and, indeed, relative number of large 
clusters in DG is greater by many orders of magnitude in comparison with the uncorrelated 
medium (see Fig. |2])p21 

At the same time Gartstein and Conwell suggested that a correlated nature of the energy 
landscape should significantly affect charge transport in disordered medium, effectively in- 
creasing the typical scale of the hop and, thus, enhancing mobility field dependence.^ The 
next major step was achieved by Dunlap, Parris, and Kenkre, who found that for charge 
carrier hopping in ID random Gaussian energy landscape^ 

/i= , P = l/ksT. (14) 

ePE J dy exp (-efiEy + /3 2 [C(0) - C{y)}) 
o 

For the case of strong dipolar disorder a^(3 3> 1 equation [TH gives 

-(a d /3) 2 + 2(<r d p)y/^0\ , (15) 



H « yu exp 



and, in general, if C(f) oc l/r p then In [i oc E p l p+l . At last, computer simulation confirmed 
that the dipolar correlation leads to the PF mobility field dependence for 3D carrier trans- 
portP^"^ (see Fig. [3} and on the basis of extended simulation an empirical relation was 
suggestecP 

H « ^exp -^dP) 2 + 0.78 [(a d /3) 3 / 2 - 2] ^eaE/a d . (16) 

Equations [15] and [16] have essentially the same dependence on E and T and only numerical 
coefficients are different. Thus, the simplified ID model could serve for a quick estimation 
of the mobility field dependence in disordered materials. 

Good qualitative description of charge transport in a correlated medium is presented in 
ref. US The principal difference between a correlated random medium and an uncorrelated 
one is that in correlated medium deepest valleys of the energetic landscape (having maximal 
escape time for a carrier at E = 0) are also the widest. For this reason carrier escape time 
for such valleys decreases drastically for E > 0. For every magnitude of E a critical size 
R C {E) exists such that valleys of this size are the most efficient for carrier capture. Simple 
calculation gives R c oc E- 1 ' 2 and R in eq [5] should be substituted for R C (E)9^ This is the 
physical reason for the PF dependence. 

An attracting feature of the DG model is absence of any significant free parameter: the 
energy scale (i.e. a) may be calculated using (P) provided p, a, and c are known. Still, 
numerical coefficients in eq[16]are certainly not constants (for example, they weakly depend 
on i?cP^), so equation [16] should be considered rather as a guide and not as a strict relation 
for fi(E,T) in any particular polar material. 

Another important feature of the DG model is that, unlike the GDM, the form of the 
mobility field dependence for moderate field does not depend on the particular kind of 
hopping rate: for eafiEjo^ < 1 mobility dependence retains its PF form (even the slope is 
the same^) for any nonpathological hopping rate. 

It is worth to mention, though, that in (partially) orient at ionally ordered dipolar matrices 
correlation properties of the distribution of C/j may be very different. In such matrices 
transport properties may have nothing in common with the usual PF picture 
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Figure 3. Mobility field dependence in the DG model for different values of 0$ (from top 
curve downward). The lowest curve is the mobility for the GDM for a/3 = 5.10. If <Jd = 0.1 
eV and a = 10A, then eaE/a d w 1 for E = 10 6 V/cm. 

CHARGE TRANSPORT IN NONPOLAR MATERIALS 

A vast collection of experimental data indicates that the model of dipolar glass cannot 
serve as a universal explanation for the PF mobility field dependence in disordered organic 
materials. If the reason for this dependence is the spatial correlation of the dipolar type fllOp . 
then it is absolutely impossible to explain the existence of the PF dependence in nonpolar 
materials, where it was routinely observed. [2i][22][3][23][2l] y[ ore c \ ose study of experimental data 
reveals, however, that in such materials the experimentally tested field range is not too 
broad (not significantly greater than one order of magnitude or even less) and, moreover, 
sometimes clearly visible deviations from straight lines can be observed when In /1 is plotted 
against E" 1 / 2 .^ In some cases these deviations lead to the upward convexity of the mobility 
curve, while in others they lead to the downward convexity. Quite formally, this behavior 
may be described by 

ln/iocE" (17) 

with n being in some cases smaller than 0.5 and in some cases greater than 0.5. This 
observation hints that, quite possibly, in weakly polar organic materials the real mobility 
field dependence is not a true PF dependence, but rather a quasi-PF dependence with 
n ^ 0.5. This quasi-PF dependence can successfully imitate true PF dependence in not so 
wide field range. According to the result of ID modelp^ the necessary condition for validity 
of relation (1171) is the algebraic behavior of correlation function C(r) oc r~ p . 
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Quadrupolar Glass Model 

One possible reason for such behavior is the contribution from randomly oriented quadrupoles. 
Quite frequently, transport molecules contain highly polar groups having significant dipole 
moments. These groups may be arranged in such a way that the total dipole moment of the 
molecule is close to zero, but the quadrupole moment of the molecule is large enough. For 
example, for the particular case of two identical groups having dipole moment p, oriented 
in opposite directions and separated by distance d, the total dipole moment of the molecule 
is exactly zero, but the quadrupole moment Q = pd is nonzero. In close analogy with the 
DG modeP we may consider the model of quadrupolar glass (QG), its simplest incarnation 
being a simple cubic lattice with sites occupied by randomly oriented quadrupoles.^ 
Calculation of the correlation function C(r) for the QG model gives^ 

0(0 "0.6^ (5)', r»a, = O (0) = ^^i-, (18) 

m m 

where c is the fraction of lattice sites, occupied by quadrupoles. If p = 3D, e = 3, a = 10 
A, b = 5 A, then for the totally filled lattice <r q ~ 0.08 eV. In ID approximation the carrier 
mobility in the QG model for the case of strong disorder (cr q /3) 2 ^> 1 has the form 

\i oc exp 

This particular field dependence can successfully imitate the PF dependence in not too 
wide field range E max /E m i n < 10. Moreover, there is another factor that can bring the 
quadrupole field dependence even closer to the PF dependence. This factor is an effect 
of dispersive transport: in dispersive regime \x oc E a W If plotted as In H vs E 3 ^, the 
dispersive component bends the straight line, making the mobility curve convex upward, 
thus pushing field dependence more close to the PF one (contribution of the dispersive 
component is evident, for example, for photocurrent transient in the data of ref. [22] even at 
high temperature). 

Detailed computer simulation of the 3D charge transport in the QG model has not been 
carried out yet, but preliminary data unambiguously support the principal result of ID model 
for the mobility field dependence ln/i oc E 3 ^ 4 '.^ Analogy with the DG model suggests that 
the functional type of the mobility dependence on E and T could be well captured by the 
ID model, while numerical coefficients should be different. Extensive comparison of the QG 
model with experimental data is still absent, though preliminary analysis of transport data 
for nonpolar polysilanes indicates that dependence In p = a + bE 3 / 4 is certainly not worse 
(in terms of statistical correlation coefficient 1Z 2 ) than the PF dependence.^ 

Los Alamos Model 

Recently a new model was suggested to explain the PF mobility field dependence in nonpolar 
polyconjugated polymers such as polyphenylene vinylens (PPVs)p2E3 According to this 
model, the major source of the energetic disorder in PPVs is almost static intramolecular 



2 3/2 



,1/2 



(19) 
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fluctuations of the torsion angle of benzene ring, resulting in the carrier energy fluctuations 
with the correlation function 

C(r) oc — exp(-ar), a = s/K, (20) 
r 

where s is the intramolecular restoring force constant and K is the intermolecular restoring 
force constant. This function in the limit case a — >■ gives an appropriate behavior needed to 
provide the PF dependence. In densely packed PPVs typical values of a should be small j^"^ 
and according to the ID model the mobility field dependence has the form 

ln/i oc -a t 2 /3 + Py/2irofeaE, o\ = v 2 /2nKa, (21) 

where a is the short-range cutoff and u is the linear electron-vibration coupling. 

Field and temperature dependence of the mobility in the LA model is different (though 
not very different) from the prediction of the QG model. Thus, a careful study of the 
dependence of \x on E and T may provide an opportunity to discriminate between these 
models. At present, there is no possibility to state which particular model better describes 
the mobility field dependence in nonpolar materials. We should note, though, that the model 
of quadrupolar glass does not relay on any particular property of the nonpolar material 
(apart from its quadrupolar nature). From this point of view it is preferable over LA model 
that, obviously, could hardly be applied for the explanation of the properties of various 
low molecular weight nonpolar organic glasses. These materials demonstrate the same PF 
field dependence as PPVs and can be perfectly well described by the QG model. High 
intramolecular flexibility is a necessary prerequisite for significant amplitude of the energetic 
disorder in the LA model, yet experimental data unambiguously demonstrate that some 
nonpolar materials consisting of very rigid molecules may still have significant energetic 
disorder (a good example of the material with planar rigid nonpolar molecules having a m 
0.09 eV is discussed in ref. I2"2l) . 

It is worth to add that the very observation that served as a first clue for the LA model, the 
significant difference between slopes of the PF dependence for two particular polyconjugated 
polymers - one is traditional MEH-PPV 311 and the other one (having very small slope of the 
mobility curve) is stiff-chain polyfluorene^ - may be reasonably well explained using the 
QG model. Indeed, the local stiffness of the main polymer chain should by itself decrease 
the qudrupolar disorder. In addition, the molecule of polyfluorene contains only carbon 
and hydrogen atoms, while the molecule of more flexible MEH-PPV contains oxygen atoms 
as well - this means that the qudrupolar moment of the monomer unit for MEH-PPV is 
greater. Hence, both factors acts in the same direction, decreasing the quadrupolar disorder 
in polyfluorene and weakening the mobility field dependence. 

Recent papers on charge transport in PPVs offer a limited opportunity to compare QG 
and Los Alamos models.^"^ These papers offer data on the temperature dependence of the 
PF coefficient &pf = j/T (see eqEj). By definition 

and for the LA model bpp oc 0. Although in the QG model coefficient bpp is not a true 
constant with respect to E, its field dependence is very weak 

6g£ oc (23) 
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and we can approximately treat it as an effective constant. Results of the fit bpp oc f3 n are 
shown in Fig. HI 

0.025 | 1 




0.000 1 1 1 1 1 1 

0.5 1 1.5 2 2.5 

n 



Figure 4. The correlation coefficient TZ 2 obtained by fitting the temperature dependence of 
6pf for different PPVs on n in bpp oc (3 n . The best fit corresponds to the minimum of 1 — TZ 2 . 
Curves for different PPVs from ref. [34] are marked A, B, C, and D; curve A' represents 
data from ref. [33] (the same PPV as for the curve A, but different thickness of the transport 
layer) . 

For most PPVs the best n is close to 1.5 and, thus, more close to 1.25 (QG model) than 
to 1 (LA model). A possible reason for the discrepancy with the QG model may be the 
specific technique used in refs. I331I341 where mobility was calculated from the current- voltage 
characteristics in the space-charge limited conduction regime assuming the PF mobility field 
dependence. Thus, the mobility field dependence was not calculated at all, only the coeffi- 
cient bpp was optimized to fit experimental current-voltage curves. The direct study of the 
temperature dependence of the corresponding true QG coefficient 



K = 7ns£ ( 24 ) 



<91n/i 

obtained from TOF experiments should be a better way to test the QG model 



Mixed Disorder: Is It Possible to Provide a Simple Formula for the Mobility Field 
Dependence? 

In a typical case the total energetic disorder in a material has several contributions of different 
nature 

s 

(say, dipolar and quadrupolar components). In such situation the correlation function (as- 
suming statistical independence of individual components) is a sum of terms having different 
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dependence on distance. The resulting mobility dependence in the ID model has the form 



l/fj, ex ^ dyexp \^-eEf3y + f3 2 ^ [C s (0) - C s (y)}j , (26) 
which cannot be written in the form of an explicit simple formula for nontrivial functions 

c s {v). 

Moreover, mobility field dependence could present a clear indication of the existence 
of several contributions to ( )25|) only when mobility is measured in an exceptionally wide 
field range. Careful analysis shows that calculation of individual dipolar and quadrupolar 
contributions to the total disorder could be performed only in the case when mobility field 
dependence is measured in the field range spanning at least three orders of magnitude.^ 
This is hardly possible in real experiments. 



TRAP-CONTROLLED CHARGE TRANSPORT 

Most organic materials contain impurities serving as traps for charge carriers (i.e. such 
impurities have energy levels laying deeper than levels of majority of transport sites). Typ- 
ically, energies of trap sites have no spatial correlation and produce an additional spatially 
uncorrelated disorder. Important question is: to what extent general properties of charge 
transport in correlated medium are insensitive to the presence of traps (later we will limit 
our consideration to the the particular case of dipolar medium)? This problem has been 
extensively studied in the ID approximation and by means of 3D computer simulation.^ 
Again, major results are the same for both approaches and may be summarized as follows: 

1) addition of traps leads to the decrease of the carrier drift mobility; 

2) PF dependence in weak fields remains unchanged without respect to trap depth A and 
concentration c; 

3) in stronger fields a linear field dependence In /i cx E emerges. 

These conclusions seem to be in striking contradiction with recent experimental observations 
for transport of holes in doubly doped polymer layersP^^ Molecules of one dopant, added 
in small concentration, and possessing significantly lower ionization potential, served as traps 
for charge carriers. In these studies it was found that for shallow traps the PF dependence 
remains mainly untouched,^"^ while for deep traps a linear field dependence was observed 
in the whole field ranged The most puzzling experimental result is the unusual dependence 
of the mobility on trap concentration JSCES 

/2 oc l/c n , (27) 

with n > 1 instead of expected dependence with n = 1 for trap-controlled transport. Again, 
this result does not agree with the theoretical dependence where n = lP^ 

Quite unexpected solution for all these puzzles was suggested in ref. [39j Surprisingly, 
the problem is not connected with any fault of the theory or inaccuracy of modern experi- 
mental technique but rather with the interpretation of the experimental data, namely with 
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Figure 5. Two methods of mobility calculation from experimentally measured photocurrent 
transient: 1) /ij = L/UE where ti is the time of intersection of asymptotes to the plateau and 
trailing edge of the transient; and 2) fii/ 2 = L/ti/ 2 E where £i/ 2 is the time for photocurrent 
to decay to the half of its plateau value. 

the particular way of the mobility calculation from the photocurrent transient temporal 
dependence. 

In TOF experiment, mobility is usually calculated by two methods (Fig. [5]). The first 
method is the method of choice for most experimental papers. Unfortunately, this particular 
method is too sensitive to the variation of the shape of the transient with variation of E (or 
c) and tends to overestimate the contribution of fast under-relaxed carriers.^ The second 
method produces mobility that is much more close to the strict definition of the mobility as 

= (v) /E where (v) is an average carrier velocity. The sensitivity of the first method is 
so significant that even in the case when good current plateau with I(t) ~ const is observed 
in the whole range of variation of the parameter (e.g., E), there is a possibility to obtain 
principally different mobility field dependences for (J>i(E) and Hi/ 2 {E) (see Fig. [6]). Note, that 
for this particular case in the whole field range the transient is more or less nondispersive and 
has a well-defined plateau (Fig. [7]). We can conclude that in some cases mobility, calculated 
by the first method, effectively masquerades a real field dependence of the true mobility. 
The same reason explains an unusual dependence of the mobility on trap concentration.^ 

The second method produces much more reliable results and should be used instead of 
the first one. This conclusion is in sharp contradiction with the common belief that field 
dependences of fii and fii/ 2 are essentially the same (experimental evidence for the difference 
may be found in literature,^ but it is usually ignored). The use of /zi/ 2 is especially important 
in situations where we have to choose among not very distinctive alternatives, such as QG 
or PF mobility field dependence in nonpolar materials. 
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Figure 6. Mobility field dependence for dipolar glass with traps for cr/3 = 3.83, A/3 = 10, 
c = 0.01, and L = 40, 000 lattice planes (equivalent to transport layer thickness of 30-40 
/im) for different methods of mobility calculation: fi v - diamonds, (iy2 - triangles, and ^ - 
squares. Here c is a fraction of sites occupied by traps. 

POOLE-FRENKEL MECHANISM REVISED: 
EFFECT OF CHARGED TRAPS 

Original PF mechanism was for a long time considered unacceptable for explanation of the 
mobility field dependence in disordered organic materials due to absence of charged traps 
(arguments provided in ref. E] suggested that trap density is certainly less than 10 13 cm -3 ). 
Yet in a recent paper the PF model was revitalized by Rackovsky and Scherj^ who argued 
that a very low density of Coulomb traps is enough to produce an essentially nondispersive 
PF charge transport (they estimated that 10 n -10 13 cm -3 should be a sufficient density for 
transport layer with the thickness of 10 /im). Because such low density can easily avoid 
detection, they suggested that the PF model still can be considered as a possible candidate 
for explanation of charge transport in disordered organic materials. 

Calculation of the PF effect in ref. 00] was carried out for the usual case of an isolated 
Coulomb trap only. Yet charged traps produce strongly spatially correlated energy landscape, 
and the hopping charge motion in such landscape usually has many features that cannot be 
captured by the approximation taking into account interaction of charge carrier with the 
isolated force center. 

Early indication that charge transport in charged medium may significantly differ from 
predictions of the model of isolated Coulomb trap was presented by Dunlap and Novikov.^ 
In recent papers^"^ this problem has been considered in more detail for a special case 
of the medium with equal concentration of randomly located positive and negative traps 
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Figure 7. Photocurrent transients for dipolar glass with traps for two values of eaE/a: 
0.021 (thin line) and 1.70 (thick line), correspondingly, with time axis rescaled to make 
values of t\i% equal for both transients. All other parameters are the same with Fig. Note 
the different shapes of transients. 



without using the isolated trap approximation.^ The result is in drastic contrast with the 
approximation of isolated Coulomb trap. For example, at E — E cr[t = 4irn(3e 3 je 1 (here n 
is trap density) a transition from mobile to immobile carriers occurs and average carrier 
velocity (v) = for E < E crit in the infinite medium, so in finite transport layers the only 
possible regime is a dispersive transport with mobility depending on thicknessP^^ 

By accident, mobility field dependence for E < E CTit retains the PF form, but its tem- 
perature dependence and, more important, concentration dependence differ drastically from 
the prediction of the isolated trap modelP^ For small concentration of traps 

ln/i « -P (c, /?) + 2 [eafi (E - £ crit ) Q] 1/2 , (28) 

3(mcr 3(lnc) cl ea 

Here the lattice version of the model is considered assuming that fraction of sites occupied 
by charges is small c< 1 but cexp h a ^> 1. 

Experimental test of the mobility dependence on the concentration of charged traps 
should be the best test of the result (|28|) . By now there are no experimental data on this 
dependence (actually, no experimental data at all on charge transport in disordered organic 
materials containing controllable concentration of static charges; quite possibly, such data 
are very difficult to obtain). Nonetheless, the mobility temperature dependence (128]) seems 
to be too strong to describe existing experimental data. For this reason we believe that the 
model of charged traps cannot be suggested as a serious candidate for explanation of charge 
transport in disordered organic materials. 
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CONCLUSIONS 



All variety of disordered organic materials could be subdivided into different classes on the 
basis of their spatial correlation properties in the distribution of energies of transport sites. 
Materials from different classes have different transport properties though differences are not 
very pronounced. In order to discriminate between different classes of the mobility field and 
temperature dependence, analysis of experimental data should be carried out with utmost 
care. In this respect an early attempt to characterize all materials in a unified manner (the 
Gaussian Disorder Model) is unjustified while we can now understand, why it was so suc- 
cessful. To some extent charge transport theory is ahead of experiment: major predictions 
(mostly concerning charge transport in nonpolar materials) still lack experimental test. For 
example, the reasonable question is: which model, quadrupolar glass or Los Alamos, is more 
suitable to describe charge transport in nonpolar PPVs? A careful and purposeful experi- 
mental study of the transport properties of nonpolar materials could resolve this problem. 
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